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Abstract 

We study the entanglement entropy Sab of a massless free scalar field on two spheres A and B 
whose radii are R\ and R2, respectively, and the distance between the centers of them is r. The 
state of the massless free scalar field is the vacuum state. We obtain the result that the mutual 
information Sa-,b = Sa + Sb — Sab is independent of the ultraviolet cutoff and proportional to the 
product of the areas of the two spheres when r»i2i,i?2, where Sa and Sb are the entanglement 
entropy on the inside region of A and B, respectively. We discuss possible connections of this result 
with the physics of black holes. 
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I. INTRODUCTION 



Entanglement entropy in the quantum field theory (QFT) was originally studied to ex- 



plain black hole entropy 



2|. Entanglement entropy is generally defined as the von Neu- 
mann entropy Sa = —TrpA In pA corresponding to the reduced density matrix pa of a 
subsystem A. When we consider the quantum field theory in (d + 1) -dimensional spacetime 
KxiV, where R and N denote the time direction and the d- dimensional spacelike manifold, 
respectively, we define the subsystem by a <i-dimensional domain A C iV at fixed time t = to- 
(So this is also called geometric entropy.) Entanglement entropy naturally arises when we 
consider the black hole because we cannot obtain the information inside the black hole. In 
fact, in the vacuum state the leading term of the entanglement entropy of A is proportional 
to the area of the boundary OA in many cases Q, This is similar to black hole entropy, 



and extensive studies have been carried out 



In this paper, we study the entanglement entropy Sab of the massless free scalar field in 
(d+l)-dimensional Minkowski spacetime on two spheres A and B whose radii are R± and R 2 
and how it depends on the distance r between the centers of the two spheres. Entanglement 
entropy of two disconnected regions has been studied, see e.g. ^-11]. We consider the case 
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that the state of the massless free scalar field is the vacuum state. We studied Sab in 
analytically. When r Ri, R2, we obtained the r dependence of Sab as 

Sab~Sa + S b , (1J 

where a is an ultraviolet cutoff length and G(Ri, R2, a) = G(R 2 , Ri,a) > 0. (Notice that we 
defined G(Ri, R 2 , a) in ([1]) as that in |l2| multiplied by ( — 1) for simplicity.) We could not 
determine the functional form of G(R\, R 2 , a). In this paper, we numerically calculate Sab 
for d = 2,3. We obtain the result that the mutual information Sa-b = Sa + Sb — Sab is 
independent of the ultraviolet cutoff length and G(Ri, R 2 , a) is proportional to the simple 
product of the surface areas of two spheres. (Note that we cannot determine the functional 
form of G(Ri, R 2 ) only by the constraints from dimensional analysis, symmetry, and behavior 
in the limit -Ri —> 0. For example, G(Ri,R 2 ) = R\R 2 + -R1-R2 is not prohibited by these 
constraints.) The mutual information is a quantity that measures the entanglement between 



two systems. (See e.g. 13[) In order to examine whether only the degrees of freedom on the 
surface of the spheres contribute to the mutual information or not, we calculate the mutual 
information So-e of two same spherical shells D and E for d = 3 and the mutual information 



Sh-i of two same rings H and I for d = 2. The internal (external) radii of the spherical shell 
and the ring are L\ (L2). The distance between the centers of the two spherical shells and 
that between the two rings are r. We obtain the result that Sd ; e and Sh-j are monotone 
decreasing function of L\. Then not only the degrees of freedom on the surface of the sphere 
but also those on the inside region contribute to the mutual information. This result is 
remarkably different from that of the entanglement entropy to which the degrees of freedom 
on the surface of the boundary contribute mainly. 



Previously, we studied Sab in UM i n or der to study an entropic contribution to the force 
between two black holes. To a distant observer, an object falling into a black hole takes an 
infinite time to reach the event horizon and the outside region is isolated from the inside 
region if we neglect the change of the mass of the black hole. Then we are probably able to 
consider the entanglement entropy of quantum fields on the outside region C of two black 
holes A and B as thermodynamic entropy, and we can see the entropic force acting on the 
two black holes from the r dependence of Sc- We consider two systems X and Y, then one 
can show Sx = Sy in general if a composite system XY is in a pure state. Then Sc = Sab 
when the state of the field on the whole space is a pure state. We will roughly estimate 
the magnitude of the entropic force between two black holes by using Sab m Minkowski 
spacetime. 



The present paper is organized as follows. There have been some computational methods 
of entanglement entropy 14 16] and the reader is urged to refer to jl7-[l9| for reviews. 



Among several others, we review in SecJTT] the method of Bombelli et al [1[ which is most 
straightforward and suitable for numerical calculations. In Sec 1111} we apply the above 
formalism to a massless free scalar field in (d + l)-dimensional Minkowski spacetime. We 
improve the computational method of Bombelli et al to reduce the computational complexity. 
In Sec JIV[ we numerically calculate the entanglement entropy Sab, the mutual information 
Sd-,e, and Sh-i in (rf+l)-dimensional Minkowski spacetime for d = 2,3. We roughly estimate 
the magnitude of the entropic force between two black holes by Sab in (3+l)-dimensional 
Minkowski spacetime. 



II. HOW TO COMPUTE ENTANGLEMENT ENTROPY 

In this section we review the computational method developed by Bombelli et al 
As a model amenable to unambiguous calculation we deal with the scalar field on M. d as a 
collection of coupled oscillators on a lattice of space points, labeled by capital Latin indices, 
the displacement at each point giving the value of the scalar field there. In this case the 
Lagrangian can be given by 

L = \G MN q M q N - \v MN q M q N , (2) 

where q M gives the displacement of the Mth oscillator and q M its generalized velocity. The 
symmetric matrix Gmn is positive definite and therefore invertible; i.e., there exists the 
inverse matrix G MN such that 

Lr LrpjV — N . [6) 

The matrix Vmn is also symmetric and positive definite. The matrices Gmn and Vmn are 
independent of q M and q M . Introducing the conjugate momentum to q M , 

Pm = GmnQ N ' i (4) 

we can write the Hamiltonian for our system as 

H = \g mn P m P n + \v MN q M q N . (5) 

Next, consider the positive definite symmetric matrix Wmn defined by 

W MA G AB W BN = V MN . (6) 

In this sense the matrix W is the "square root" of V in the scalar product with G. 

Now consider a region Q in M. d . The oscillators in this region will be specified by Greek 
letters, and those in the complement of Q, Q c , will be specified by lowercase Latin letters. 
We will use the following notation 

(w ab W a0 \ ( A B\ 4R (w ah W als \ ( D E\ 

W AB = aP = W AB =\ = , (7) 

\W ab W aP ) \B T C J \W ab W«P J \E T F) 

where W AB is the inverse matrix of Wab {W ab is not obtained by raising indices with G AB ). 
So we have 

/ AD + BE T AE + BF 





B T D + CE T B T E + CF 



(8) 
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If the information on the displacement of the oscillators in Q is considered as unavailable, 
we can obtain a reduced density matrix p red for Q c , integrating out over q a e R for each of 
the oscillators in the region Q, and then we have 

Predtf, q' b ) = I I! tfPtf' <A <A ^ ( 9 ) 

where p is a density matrix of the total system. 

We can obtain the density matrix for the ground state by the standard method, and it 
is a Gaussian density matrix. Then, p Te d is obtained by a Gaussian integral, and it is also a 
Gaussian density matrix. The entanglement entropy S = —trp red In p red is given by [l| 

S = ^2f(K), (10) 

n 

/(A) = ln^A 1 ' 2 ) + (1 + A) 1/2 ln[(l + A" 1 ) 1 / 2 + A" 1 / 2 ], (11) 
where A n are the eigenvalues of the matrix 

In the last equality we have used dSJ. The last expression in ( !T2|) is useful for numerical 
calculations when f2 c is smaller than Q, because the indices of A and D take over only the 
space points on Q c and the matrix sizes of A and D are smaller than those of B and E. It 
can be shown that all of A n are non-negative as follows. From (|SJ) we have 

AA = —AEB T = BFB T . (13) 

It is easy to show that A, C, D and F are positive definite matrices when W and W~ Y are 
positive definite matrices. Then AA is a positive semidefinite matrix as can be seen from 
( IT3l) . So all eigenvalues of A are non-negative. After all, we can obtain the entanglement 
entropy by solving the eigenvalue problem of A. 



III. LATTICE FORMULATION 



We apply the above formalism to a massless free scalar field in (d + l)-dimensional 
Minkowski spacetime. The Lagrangian is given by 



d* x \[4> 2 -<y4f\. (14) 



As an ultraviolet regulator, we replace the continuous (^-dimensional space coordinates x by 
a lattice of discrete points with spacing a. As an infrared cutoff, we allow the individual 
components of n = x/a to assume only a finite number N of independent values —N/2 < 
< N/2. The Greek indices denoting vector quantities run from one to d. Outside this 
range we assume the lattice is periodic. The dimensionless Hamiltonian H = aH is given 
by 

1 1 ^ 2 2 1 1 

H = aH = J^l^n + 2 J2^+ s ^ " ^ + ^T^] = H 2 71 "" + S 2^ mVmn ^ (15) 

n ^=1 n m,n 

where (f) n and 7r n are dimensionless and Hermitian, and obey the canonical commutation 
relations 

[<f>n, TTm] = «<W (16) 

In Eq.( !T5|) we insert a mass term in order to remove a zero eigenvalue of V mn ; if V mn should 
have the zero eigenvalue, W~ l in (J7J) would not exist. Later we will take N to infinity. In 
this limit we can neglect the zero eigenvalue of V mn and will take am to zero. Taking N 
to infinity is important in order to calculate the entanglement entropy Sab of two spheres. 
The entanglement entropy of two spheres is more sensitive to the value of N than that of 
one sphere. (In fact, we numerically calculated Sab for finite iV with antiperiodic boundary 
conditions without the mass term. Sab depends on N when the distance r between two 
spheres is close to N/2, and we could not obtain the clear r dependence of Sab-) 
From ffTol) we obtain (see e.g. |20|) 

_ d 9 —u 

W mn = N~ d Yy™ 2 + 2 J> - cos =^)} 1 / 2 e 2mk ^ n - m ^ N , (17) 

k fi=l 

W~l = N~ d J> 2 m 2 + 2 £(1 - cos H^-V^™)/* (18) 

k /i=l 

where the index k also carries d integer valued components, each in the range of —N/2 < 
k^ < N/2. We take N to infinity and change the momentum sum into an integral with the 
replacements = 2-nk^/N and N~ d — > 7^3, and then we have 

id. 



W mn = f A±- e ^- m \a 2 m 2 +2^(1- cos g M )]5, (19) 

w ™ = £ ^ eiq(n ~ m) ^ m2 + 2 X> - cos ^)]^- (2°) 
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In ( TT9|) and ( 120]) the integrals converge when am — > 0, so we can take am to zero, 

W ™ = {(l^© 1 ^ 08 ^^ (21) 

= £ ^/ q{n ' m ^YM ~ «*<&.)]^- (22) 

From (12T|) and ( 122|) we can compute W 7 ™™ and W 7 "* numerically. Then we can compute 
the entanglement entropy from ffTU]) . f|TT]) and f|T2|) . The integrands in fl2Tl) and (1221) highly 
oscillate when \\n — m\\ ^> 1, and the numerical integrals converge very slowly We can 
obtain approximate expressions of W mn and by hand when \\n — m|| > 1, so we will 
use them when \\n — m\\ ^> 1 in order to reduce the computational complexity of W mn and 
W~^. To evaluate W mn and when \\n — m\\ ^> 1, we define r = a(n — m) and take 
\\n — m\\ to infinity keeping r fixed. We change the variable as p = q/a, and then we have 



/_! ra^PS 1 - C08 *>^ - «*" /_°° tSj^^DWII*. (23) 

a M — 1 



We can perform the integral in (1231) analytically when ||r||/a — > oo (see Appendix A of 12] 
), and then we obtain 



where 



(d - 1)!! 

for even d > 2, 



*=< ( tA" (25) 

-2 / " for odd d > 3. 

(27r)( d+1 )/ 2 



We can evaluate W m ^ when ||n — m|| 3> 1 in the same way (see Appendix A of |l2j ), and 
then we obtain 

^ °" £ ^>'-' m ^ - = j^p. (26) 

—. — for even d > 2, 

2 (2^mv^ for odd d " 3 ' 

where 0!! = (-1)!! = 1. 



where 



IV. NUMERICAL CALCULATIONS 



We calculate numerically the entanglement entropy Sab of two spheres A and B whose 
radii are R\ and R 2 , and the distance between the centers of them is r for d = 3. 

We put the centers of the spheres on a lattice. We define the sphere whose radius is 
R as a set of points which are at distances of R or less from the center of the sphere. In 
order to reduce the computational complexity of W mn and W~^, we use the approximate 
expressions (T2"4"j) and (126]) when ||n — m\\ > 10, and we use the numerical integrals of exact 
expressions (12T!) and (|22|) when \\n — m\\ < 10. When \\n — m\\ = 10, the differences between 
the numerical integrals of the exact expressions and the approximate expressions are less 
than 4% for W mn and less than 1% for W~^. We perform matrix operations and calculate 
the eigenvalues X n of the matrix A in f[T2T) with Mathematica 8. The number of columns 
and rows of A is the number of points in the region of which we calculate entanglement 
entropy. 

We show the computed values of S(R) which is the entanglement entropy of one sphere 
as a fumction of R 2 /a 2 in FigJT], where a is a lattice spacing. The points are fitted by a 
straight line: 

S = 0.37R 2 /a 2 . (28) 

This result agrees with the result in except for the coefficient. (The coefficient in is 
0.30. This difference necessarily arises from the difference of regularization methods. In 
the author use the polar coordinate system and replace the continuous radial coordinate by 
a lattice.) 

We show the computed values of Sab(i", Ri, R2) which is the entanglement entropy of two 
spheres as a function of r/a for R 1 /a = R 2 /a = 6, 7 in FigO As can be seen, Sab reaches 
its maximum value Sa + Sb when r — > 00. In order to clarify the behavior of Sab as a 
function of r, we show the computed values of (Sa + Sb — Sab) -1 ^^, Ri, R2) as a function 
of r/a for Ri/a = R2/C1 = 6, 7 in Figj3l The straight lines in Fig 13] are fitted by the data 
between r/a = R\/ 'a + Rij 'a + 24 and r/a = R\ /a + Ri/a + 84. In these regions the points 
are beautifully fitted by the straight lines. Then, when r 3> R\, R2, we obtain 

- S A ; B = S AB (r, R,, R 2 ) - S A (Rx) - S B (R 2 ) » _ g ( i ^ i ^) > ( 29 ) 

where G(i?i,i?2) is defined in (129]) and G(Ri,R2) = G(R2,Ri) > 0. Sa-b is the mutual 



S(R) 




R 2 /a 2 



FIG. 1. The entanglement entropy S(R) of one sphere whose radius is R as a fumction of R 2 /a 2 . 
The line is the best linear fit. 

information of A and B. From FigJ3]the approximate expression (1291) is precise for relatively 
small r. (When Ri = R 2 = R, for r <J 3R (|2"9"1) is precise from FigJ2J) We can obtain 
G(Ri, i?2)/a 4 from slopes of graphs of (Sa + Sb — Sab)' 1 ^^, Ri, R2), and then we show 
the computed values of G(R\, R 2 )/a A as a function of R\/a 2 for Ri/a = 4, 4. 5,..., 7 in 
FigJH From FigHJwe can see that G(Ri, R 2 )/a A is proportional to R\. Because G(Ri, R 2 ) = 
G(R 2 , Ri), we obtain G(Ri,R 2 ) = gR 2 R 2 , where g is a dimensionless constant. We can 
obtain the values of gR\ from slopes of graphs of G(Ri, R 2 ) as a function of R\ . To obtain 
the precise value of g, we show the computed values of gR\/a 2 as a function of R\/a 2 in 
Figj5]and obtain g = 0.26 from the slope of the line which is the best linear fit in Figj5j 
Finally, when r 3> Ri, R 2 , we obtain 

0.26R 2 R 2 



Sa-,b = Sab^-, Ri, R2) — Sa(Ri) — Sb{Ro) 



(30) 



When r Ri R 1 ,R 2 , from Fig|2J Sab rapidly decreases when r decreases. (Note that we 
cannot determine the functional form of G(R\, R 2 ) only by the constraints from dimensional 
analysis, symmetry, and behavior in the limit Ri — > 0. For example, G(Ri,R 2 ) = R\R 2 + 
RiR\ is not prohibited by these constraints.) 

For d = 2, we compute Sab in the same way. We show only the computed values of 
(Sa + Sb — Sab^ 1 ^ 2 ^', Ri, R2) as a function of r/a for R 1 /a = R 2 /a = 15,16 in FigEJ 
The straight lines in FigJ6] are fitted by the data between r/a = Ri/a + R 2 /a + 101 and 
r/a — Ri I a + R 2 /a + 201. In these regions the points are beautifully fitted by the straight 
lines. We cross-checked our numerical procedure with the data of related calculations in 



Figure 1 in 



(In the figure, the authors show the mutual information of two discs for 
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FIG. 2. Sab — Sa~ Sb as a function of r/a for R\/a = R 2 /a = 6, 7, where Sab is the entanglement 
entropy of two spheres A and B whose radii are R\ and R 2 . The distance between the centers of 
the two spheres is r. 




FIG. 3. (Sa + Sb — Sab)~ 1 ^ A as a function of r/a for R\/a = R 2 /a = 6,7, where Sab is 
the entanglement entropy of two spheres A and B whose radii are R\ and i?2- The distance 
between the centers of the two spheres is r. The straight lines are fitted by the data between 
r/a = Ri/a + R 2 /a + 24 and r/a = R 1 /a + R 2 /a + 84. For r > 3R{= R x = R 2 ) the lines are 
beautifully fitted and the approximate expressions f|29[) and (|30p are precise. 

Ri = R 2 = R and r = 3i?. Our results were very close to theirs.) Finally, when r ^> Ri, R 2 , 
we obtain 

- S A; b = S AB (r, R x , R 2 ) - S A (Ri) - S B (R 2 ) w _ °- 37 ^ R \ (31) 

In order to examine whether only the degrees of freedom on the surface of the spheres 
contribute to the mutual information or not, we calculate the mutual information Sd-,e of two 
same spherical shells D and E for d = 3 and the mutual information Sh-i of two same rings 
H and / for d = 2. The internal (external) radii of the spherical shell and the ring are L x 
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FIG. 4. G(Ri,R 2 )/a A in ([29]) as a function of ^/a 2 for R x /a = 4,4.5, ... ,7. The lines are the 
best linear fit. 

gRiV 




FIG. 5. gR 2 /a 2 as a function of R 2 /a 2 , where g is defined as G(R\,R 2 ) = gR\R\- The line is the 
best linear fit. 

(L 2 ). The distance between the centers of the two spherical shells and that between the two 
rings are r. When r 3> L 2 , we obtain Sd-,e ~ G ss (Li, L 2 )/r 4 and SV;i ~ G r (Li, L 2 )/r 2 . We 
show (G„(Li, L 2 )) X I 2 IL 2 for L 2 = 10a as a function of LJL 2 in FigHand (G r (Li, L 2 )) 1 / 2 /L 2 
for L 2 = 22a as a function of L x jL 2 in FigE The curve in Figdis 0.50(1 - (L 1 /L 2 ) 3 ) 2 / 3 
and the curve in Figj8]is 0.56(1 — (Li/L 2 ) 2 ) 1 ^ 2 . We show these curves for comparison with 
the data. From FigHand FigJH (G S8 (Li, L-2)) 1 / 2 / L 2 2 and (G r (Li, L 2 )) 1/2 / 1 L 2 are monotone 
decreasing function of Li/L 2 , and (G ss (Li, L 2 )) 1 / 2 is not proportional to the 2/3 power of 
the volume of the spherical shell, and (G r (L 1; L 2 )) 1 / 2 is not proportional to the 1/2 power 
of the area of the ring. Then not only the degrees of freedom on the surface of the sphere 
but also those on the inside region contribute to the mutual information, and the degrees of 
freedom on the inside region does not contribute uniformly to the mutual information. 
We roughly estimate the magnitude of the entropic force between two black holes by using 
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FIG. 6. (Sa + Sb — Sab) _1//2 as a function of r/a for R\/a = R 2 /a = 15,16, where Sab 
is the entanglement entropy of two discs A and B whose radii are R\ and R2. The distance 
between the centers of the two discs is r. The straight lines are fitted by the data between 
r/a = Ri/a + R 2 /a + 101 and r/a = R\/a + R 2 /a + 201. For r > 4R(= R 1 = R 2 ) the lines are 
beautifully fitted and the approximate expression (|31 1) is precise. 

(G ss (L,,L 2 )) 1/2 /(L 2 ) 2 




FIG. 7. (G ss (L 1 ,L 2 )) 1/2 /L 2 2 for L 2 = 10a function of L\/L 2 . G ss (Li, L 2 ) is defined as 

S D ;E ~ G ss (L!,L 2 )/r 4 when r » L l5 L 2 . The curve is 0.50(1 - (Z^/La) 3 ) 2 / 3 . {G SS (L U L 2 )Yl 2 /L\ 
is monotone decreasing function of L\/L 2 and not proportional to (1 — (Li/L 2 ) 3 ) 2 ^ 3 . 



Sab in Minkowski spacetime. We consider two black holes (A and B) which have the same 
radius R\ = R 2 = R and the distance between which is r. For simplicity, we consider the 
case that the state of the field on the whole space is a pure state. Generally, if a composite 



system XY is in a pure state, then Sx = Sy [131 ] . Then the entanglement entropy of the 
outside region of two black holes is equal to that of the inside regions of two black holes. 
We define the entropic force of the field on the outside region which acts on one black hole 
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FIG. 8. (G r (Li, L2)) X I 2 IL\ for L 2 = 22a as a function of Li/L 2 . G r (Li,L 2 ) is defined as 
S H -i ~ G r (Li,L 2 )/r 2 when r » Li,L 2 . The curve is 0.56(1 - (Li/L 2 ) 2 ) 1/2 . (G r (Li, Z^)) 1 / 2 /^ is 
monotone decreasing function of L\/L 2 and not proportional to (1 — (L\/ L 2 ) 2 ) l l 2 . 



in the direction of increasing r as F e f. F e f is given by 

„ rr,dS AB 



ef 



dr 



(32) 



where T is the temperature of the field of the outside region. To estimate F e f, we set Sab 
to that in Minkowski spacetime and T to the Hawking temperature T = (SttGnM)' 1 = 
(4irR)~ l . In this approximation the entropic force is repulsion force because S A b increases 
when r increases. dS A B / dr is independent of the ultraviolet cutoff, and then we obtain 



F ef = ~S' A . >B (r/R) 



AnR 2 



S' A . B (r/R), 



(33) 



where S A] b = S A (R) + Sb(R) — Sab^, R) and S' A . B = dSA-,B/d(r/R). (Sa-b is independent 
of the ultraviolet cutoff and a function of r/R.) Then the ratio of the entropic force to the 
force of gravity (F g = -^f- = -jS^) is 





Fef 




S'a;b{t/R) 




F 9 




R 2 /r 2 


where lp is the Planck length I 


p = 


(G N h/c 3 ) 1/2 . 


When r ^> 



(34) 



and then we obtain 



, lp\ 2 R 3 
0.33 1^1 — 



(35) 



When r « R, we show the computed values of (j^j ~p~ as a function of r/R for i?/a = 10 



in Figj9j (Although from fl34]) 



r 


Fef 




F 9 



is a function oi r/R and is independent of the 



choice of the value of R/a, the computed values of 



r 


Fef 




F 3 



slightly depend on the choice 
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(RJl P f*(FJF g ) 
5+ 



2.1 



2.2 2.3 2.4 2.5 



2.6 2.7 



r/R 



to 






F 9 



as a function of r/R for 



i2/o = 10. 



of the value of R/a because the spheres on the lattice are distored. When R/a is large, the 
spheres on the lattice are similar to the real spheres and this R/a dependence is small.) From 
(|35|) and Fig J9] the entropic force is much smaller than the force of gravity when i? 3> lp and 
comparable to the force of gravity when R m lp. 

V. CONCLUSIONS AND DISCUSSION 



We calculated numerically the entanglement entropy Sab of two spheres and obtained 
the approximate expression ( 1301) . From Figj3j ( 1301) is precise for relatively small r. (When 
R\ = i?2 = R, for r ^ 3i? ( 130]) is precise from Figj3j) We showed that the mutual information 
Sa-b of ^4 and B is independent of the ultraviolet cutoff for d — 2, 3 though 5^ and 
depends on the ultraviolet cutoff. The mutual information Sa,b measures the entanglement 
between A and B and Sa measures the entanglement between A and A c where A c is the 
complementary of A. Then our results mean that the ultraviolet divergence of entanglement 
entropy in QFT is caused by the entanglement between points which are infinitely close to 
each other and the entanglement between regions which are finitely separate from each other 
is finite. And we showed that Sa-b is the simple product of a function of R± and that of 
R 2 for d = 2, 3. These properties of Sa-,b for d = 2, 3 are most likely the same as those for 
d > 4. Then, from (DO), for d > 4 when r>i?i,i?2 we assume 



S 



A\B 



where ga > is a dimensionless constant. 



gdRj 1 R I 2 1 

r 2d-2 
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(36) 



In order to examine whether only the degrees of freedom on the surface of the spheres 
contribute to the mutual information or not, we calculate the mutual information Sd-e of 
two same spherical shells D and E for d = 3 and the mutual information Sh-,i of two same 
rings H and / for d = 2. We obtained the result that not only the degrees of freedom 
on the surface of the sphere but also those on the inside region contribute to the mutual 
information, and the degrees of freedom on the inside region does not contribute uniformly to 
the mutual information. Because Sd-,e and Sh-j measure the entanglement between regions 
which are finitely separate from each other, it is natural that the inside region contribute to 
the mutual information. The result that the inside region does not contribute uniformly to 
the mutual information means that the mutual information is not the product of the simple 
sum of the contribution from each volume elements. These results are different from that of 
the entanglement entropy to which the degrees of freedom on the surface of the boundary 
contribute mainly and uniformly. So the mutual information of two disconnected regions 
is not universally proportional to the product of the surface areas of the regions. Because 
a sphere has only one dimensionful parameter, the mutual information of two spheres is 
proportional to the product of the surface areas. For example, the mutual information of 
two rectangular solids is most likely not proportional to the product of the surface areas 
because a rectangular solid has three dimensionful parameters. 

Our numerical method has three properties. First, we take the volume of the whole 
space to infinity, i.e. N — > oo in ffT7|) and ffT8]) . Second, the computational complexity of 
our method depends only on the number of points on the regions of which we compute the 
entanglement entropy and does not depend on the distance between the separated regions. 
The computational complexity of conventional methods increases when the distance between 
the separated regions increases. This is because the numerical integrals of W mn in (T2~Tj) and 
in fl22|) converge very slowly when \\n — m\\ ^> 1. In order to reduce the computational 
complexity of W mn and W~^, we use the approximate expressions ()2l]) and (126|) when \\n — 
m\\ > 10. Third, we can compute the entanglement entropy of general shaped regions 
by our method because we do not use any symmetry of the regions of which we compute 
the entanglement entropy in our method. For example, we can compute the entanglement 
entropy of more than two separated regions. The first and the second properties enable us 
to obtain the r dependence of Sab- And the third property enable us to compute Sab f° r 
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We estimated roughly the magnitude of the entropic force between two black holes. From 
fl35|) and Figj9]the entropic force is comparable to the force of gravity when R w lp. This 
rough estimate suggests that the entropic force is important for Planck scale black holes. 
(Of course, this result would be changed if the effect of quantum gravity would be taken 
into account when R f» l P .) 

Next, we discuss the microscopic origin of the entropic force. As we see from (1331) the 
entropic force is proportional to the r derivative of the mutual information Sa-b- So the 
origin of the entropic force is the entanglement between inside regions of two black holes. 
Due to the entanglement between inside regions of two black holes, the density matrix of 
the scalar field on the outside region changes when r changes. Then the force acts on black 
holes along the direction in which Sab increases. 

Finally we mention the validity of this estimate. When r 3> R, it is shown that Sa-b in 
the black holes case can be expected to be similar to that in the Minkowski spacetime case 
except for the coefficient because almost all regions between two black holes is similar to 
Minkowski spacetime [ijj]. So, the rough estimate corresponds to the contribution to F e f in 
fl32|) from Sa-,b- However, in the black holes case Sa and Sb depend on r and contribute to 
F e f. These contribution from Sa and Sb has been discussed in fl^j ]. When r w R, Sa-,b in 
the black holes case is probably different from that in the Minkowski spacetime case because 
the region between two black holes is very different from Minkowski spacetime. However, 
even when r m R, Sa-,b is most likely independent of the ultraviolet cutoff as that in the 
Minkowski spacetime case. 
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